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MULTI-CHANNEL FREQUENCY DOMAIN EQUALIZER 
FOR RADIO FREQUENCY INTERFERENCE CANCELLATION 

Background of the Invention 

1. Field of the Invention 

This invention relates generally to radio frequency interference impairments 
associated with digital subscriber lines, and more particularly to a method for exploiting 
deterministic spreading of the fast Fourier transform (FFT) at the receiver to cancel 
correlated noise sources such as RFI. 

2. Description of the Prior Art 

ITU Telecommunications Standardization Sector, G.test.bis., Study Group 4, 
Question 15, 1999, has proposed various radio frequency interference (RFI) impairments 
for incorporation into digital subscriber line (DSL) performance tests. For asymmetric 
DSL (ADSL) systems such as described in ITU Telecommunications Standardization 
Sector, G.992.1: Asymmetrical Digital Subscriber Line (ADSL) Transceivers., Study 
Group 4, Question 15, 2000, RFI is assumed to be caused by narrowband amplitude 
modulated (AM) radio stations, which overlap the downstream ADSL band from 540- 
1 104 kHz. Additional RFI sources include, but are not limited to, aliases from AM 
stations operating above 1104 kHz, amateur (HAM) transmitters, televisions, computer 
equipment, and various other (relatively) low frequency signals (e.g., aeronautical 
radiolocation signals). 

At the ADSL transmitter, data is assembled into subchannels, transformed by an 
inverse fast Fourier transform (IFFT), and a cyclic prefix is added. After propagating 
through a communication channel, the opposite operations are performed by the receiver. 
The cyclic prefix is removed, and the received data is transformed by the fast Fourier 
transform (FFT). If the communication channel is shorter than the cyclic prefix and the 
noise is uncorrected on the individual subchannels, then the transmitted data can be 
recovered by a one-tap equalizer (the frequency-domain equalizer or FEQ), such as 
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described by T. Starr, J. Cioffi and P. Silverman, Understanding Digital Subscriber Line 
Technology, Upper Saddle River, NJ: Prentice Hall, 1999. 

The noise that results from RFI at the receiver is deterministically spread by the 
FFT to neighboring subchannels. As a result, instead of one (or a few) subchannels being 
effected by the narrowband impairment, many subchannels are effected. The resulting 
loss of signal-to-noise ratio (SNR) can unacceptably limit the capacity of the 
communication channel 

In view of the foregoing, a need exists for a scheme to exploit the deterministic 
spreading of the FFT at the receiver to cancel correlated noise sources such as RFI in 
order to increase the effective data communication rate between subscribers using shared 
media. 
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Summary of the Invention 

The present invention is directed to a method of estimating the multichannel FEQ 
in block and iterative fashions. Specifically, algorithms are developed to estimate the 
multichannel FEQ. The estimated multichannel FEQ provides improvements in SNR 
over that achievable using standard single channel FEQ, 

According to one embodiment, a multichannel FEQ is formulated when there is 
correlation in the noise matrix (e.g., from a correlated noise source such as RFI). 
Equalizers which have nonzero elements in addition to the main diagonal are considered. 
Each row of the multichannel FEQ is recovered via choosing the nonzero columns in 
addition to the column directly correlated with the input through the channel. Columns 
corresponding to subchannels having strong correlations with the column directly 
correlated with the input through the channel are most preferably chosen (generally 
neighboring subchannels, or subchannels corresponding to strong RFI sources). As the 
RFI and signal power become closer, not transmitting data on a subchannel allows that 
subchannel to be used more effectively with other subchannels in the multichannel FEQ. 

As used herein, FEQ means a complex multiply to adjust for the gain and phase of 
the effective channel response in each sub-symbol. The FEQ does not enhance the SNR; 
it only scales and rotates each sub-symbol so that a common decoder can be used. 

As used herein the following notation will have the following meanings. Italics 
will be used for scalars, and boldface for vectors and matrices. A superscript ~ l will refer 
to the pseudoinverse, * to conjugation, to transpose, and to conjugate transpose. The 
operator diag() forms a diagonal matrix with the arguments being the elements of the 
main diagonal, wherein E[] is the expected value of its argument, and Amax() is the 
maximum eigenvalue of its argument. The term a:b will be used to denote a contiguous 
subset of a vector or dimension of a matrix, and a ":" by itself refers to all of the elements 
along that dimension. For example, if 

c = [c(l)...c(L)] r , 

then, 
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c(:) = c, and 

c(a:b)-[c(a)...c(b)] r 

To select a (possibly) noncontiguous subset of a vector or dimension of a matrix, a vector 
of the desired indices will be used as the argument of that dimension. For example, if 

k = [£(l)...^M)f 
is a vector of integers, then 

c(k)~[c(k(l))...c(k(M))f. 
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Brief Description of the Drawings 

Other aspects and features of the present invention, and many of the attendant 
advantages of the present invention, will be readily appreciated as the same become 
better understood by reference to the following detailed description when considered in 
connection with the accompanying drawing wherein: 

Figure 1 is a graph illustrating a simulated subchannel SNR for additive white 
Gaussian noise (AWGN) with no RFI, AWGN with RFI using single channel FEQ, and 
AWGN with RFI using multichannel FEQ in accordance with one embodiment of the 
present invention; 

Figure 2 is a system diagram depicting a multichannel FEQ integrated with an 
ADSL communication system according to one embodiment of the present invention; and 

Figure 3 is a flow diagram illustrating a method of estimating multichannel FEQ 
coefficients for the n th subchannel, g(w), according to one embodiment of the present 
invention. 

While the above-identified drawing figures exemplify characteristics associated 
with particular embodiments, other embodiments of the present invention are also 
contemplated, as noted in the discussion. In all cases, this disclosure presents illustrated 
embodiments of the present invention by way of representation and not limitation. 
Numerous other modifications and embodiments can be devised by those skilled in the 
art which fall within the scope and spirit of the principles of this invention. 
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Detailed Description of the Preferred Embodiments 

The present invention is best understood by considering a system with n = 1, . . ., 
N subchannels which transmits data at time t = 1, . . T. If subchannels are arranged 
along the rows and times are arranged along the columns, then a Nx T transmitted data 
matrix X can be written as 



X = 



X(l,l) X(IJ) 



X(N,1) - X(NJ) 



(1) 



Y, a Nx T matrix of received symbols, and V, a A r x T noise matrix, are defined in an 
analogous manner. 

In order to implement one embodiment of the present invention, the present 
inventor designed aiVxJV equalizer, G, such that 



G(n, :)= argmin£[lX(^/)-G(« J :)Y(:,OfJ 



(2) 



G(k,:) 



If the channel matrix H = diag (H(l) . . . H(Nj)> and if the received data at time t can be 
written as 



Y(:,0 = HX(:,0 + V(:,f), 



(3) 



where the elements of V(:, 0 are uncorrected, then the zero forcing equalizer G can be 
recovered as 



G = diag 



f i i ^ 



H(l)'-' H(N)J 



(4) 



The minimum mean-square error solution also has a diagonal structure: 
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G = diag 



mr 



\HiX)\ 



2 Px{l)' 

+ - 



Pv(l) 



H(N)* 

\H(N) \ 2 +£^ll 
1 A P V (N)) 



(5) 



where Px(k) is the power spectral density (PSD) of the input, and Pv(k) is the PSD of 
the noise on the k* subchannel. These arguments are seen to be the standard one-tap 
FEQ's such as described by T. Starr, J. Cioffi and P. Silverman, Understanding Digital 
Subscriber Line Technology, Upper Saddle River, NJ: Prentice Hall, 1999, referenced 
herein before. 

The present inventor recognized, however, when there is correlation in the noise 
matrix V (e.g., from a correlated noise source such as RFI), it seems reasonable to 
consider equalizers G that have nonzero elements in addition to the main diagonal. One 
embodiment can be exemplified by considering the n th subchannel and letting k(n) be the 
indices of the desired nonzero columns of the n th row of G. Most preferably, G(n, k(n)) 
can then be found such that 



X(«,0 = Y(k(«),O r G(«,k(n)) 7 



(6) 



Letting g(w) = G(«, k(«)) r , .collecting data from time t = 1, . . ., T, and making the obvious 
substitutions then yields: 



" X(n,V)~ 




"Y(k(n),l) r " 


X(n,T) 




_Y(k(n),J)^ 



g(«): 



(7) 



where 



x(«) = Y(n)g(«) 
Therefore, the n th row of the equalizer G can be recovered as 



(8) 
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g(n) = Y(«)-x(«) 



(9) 



This procedure is then repeated for the remaining subchannels to form the complete 
multichannel FEQ: 



G = 



g(l) 7 " ' 



(10) 



The present invention is not so limited however, and it shall be understood that 
other solutions, including but not limited to, iterative solutions are also possible using 
least mean-squares (LMS), recursive least-squares (RLS), and their related variants such 
as discussed by G.-O. Glentis, K. Berberidis and S. Theodoridis, "Efficient least squares 
adaptive algorithms for FIR transversal filtering," IEEE Signal Processing Magazine, 
vol. 16, no. 4, pp. 13-41, 1999. The n th row, for example, of the multichannel FEQ can 
be updated as: 



e(i) = X(n,t)-Y(k(n)jfg(n) 



(11) 



g(«) = g(») + //(fM0Y(k(n),0* 



(12) 



for t= 1, . .., T, where pi(t) controls the adaptation. Popular choices include: 



Standard LMS: 



(13) 



where 
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0< a <_!_ (14) 

X max(R) 



and 

R = £[Y(k(n),0Y(k(«),0"] (15) 

• Normalized LMS (NLMS): 

M(t) = p + YiHnltfYikinlt)' ^ 
where a e (0,2) and 0 < j3 . 

• Power Normalized LMS (PNLMS): 



»v=ikY (17) 



2 2 

where cr 2 (r) = ccr 2 (f-l) + |e(0| , ce (0,1), and 0<a< — . 



Most preferably, when choosing the nonzero columns of the n th row of G ? other 
than the n h column (which should be selected as this is the column directly correlated 
with the input through the channel), columns corresponding to subchannels with strong 
correlations with the n th subchannel should be chosen. These are often neighboring 
subchannels, or subchannels corresponding to strong RFI sources. 

Subchannels with strong RFI can be found via any well-known peak picking 
algorithm applied to the noise vector V(:, f). It can be appreciated that when RFI 
dominates the received signal on a subchannel, then adding a weighed version of the 
received subchannel with RFI to the desired subchannel n is (almost) equivalent to the 
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same operation when no data was transmitted on that subchannel. As the RFI and signal 
power become closer, however, not transmitting data on a subchannel allows that 
subchannel to be used more effectively with other subchannels in the multichannel FEQ. 

A good model for the RFI (or other source of subchannel noise correlation) could 
potentially be exploited to improve the performance such as discussed by B. Wiese and J. 
Bingham, "Digital radio frequency cancellation for DMT VDSL," ANSI T1EL4/97-460, 
Dec. 1997. A positive aspect of the present method however, is that its adaptive nature 
allows it to be somewhat decoupled from the exact RFI model. Further, the single 
channel FEQ or an existing RFI model could be used to initialize the present adaptive 
process. 

Table 1 below highlights the channel capacity in Mbps for three different 
scenarios including AWGN with no RFI using a single channel FEQ, AWGN with RFI 
using a single channel FEQ, and AWGN with RFI using the present multichannel FEQ. 
The results were obtained by simulating a downstream ADSL channel with RFI added at 
the receiver. Three subchannels were used for the design of the multichannel FEQ for the 
n th subchannel: the n th subchannel, and the two other subchannels with the strongest noise 
power. The multichannel FEQ coefficients were solved for using the NLMS algorithm. 
The channel was modeled as 6000 feet of 26 AWG twisted pair copper wire with -140 
dBm/Hz additive white Gaussian noise (AWGN) and RFI model 1 impairments such as 
discussed by ITU Telecommunications Standardization Sector, G.testbis, Study Group 4, 
Question 15, 2000, referenced herein before. The capacity (Q in bits/s (bps) was 
calculated from: 

C = /olog 2 Q + £NR), 
where f 0 is 4000, and SNR is the gap adjusted signal-to-noise ratio after the FEQ. 



Patent Application of Arthur J. Redfern 



Page 10 of 18 



Docket No,: TI-3 1367 



Table 1: Capacity calculations 



Noise 


FEQ 


Capacity (Mbps) 


AWGN 


Single channel 


10.951 


AWGN+RFI 


Single channel 


5.794 


AWGN+RFI 


Multichannel 


7.681 



Figure 1 is a graph illustrating a simulated subchannel SNR for 1) additive white 
Gaussian noise (AWGN) with no RFI using a single channel FEQ, 2) AWGN with RFI 
using a single channel FEQ, and 3) AWGN with RFI using the present multichannel FEQ 
in accordance with one embodiment of the present invention. The results suggest that 
significant SNR improvements (and corresponding capacity gains) are possible with the 
present multichannel FEQ. 

Figure 2 is a simplified block diagram depicting an ADSL communication system 
10 using a multichannel FEQ 12 according to one embodiment of the present invention. 
At the transmitter, input data 14 to be transmitted is segmented (serial-to-parallel (S/P) 
converted), as depicted in element 16, into blocks of length N. The n th element is referred 
to as the n th subchannel. The data is transformed (e.g., in the case of ADSL, by an IFFT 
as depicted in element 18, and the addition of a cyclic prefix as depicted in element 20) to 
length N + P. The transformed data is next unsegmented, as depicted by element 22 
(parallel-to-serial (P/S) converted). The data is then sent through the channel as shown 
by element 24 (perhaps after additional filtering, mixing, etc. not shown). This channel 
data is then received by the ADSL receiver, perhaps after filtering, mixing, etc., and 
generally includes time domain equalization (TEQ) shown by element 26, that is 
conventionally implemented as an FIR filter. Noise 28, such as RFI discussed herein 
before, may be introduced into the signal before it is received by the ADSL receiver. The 
noise that results from RFI at the receiver is deterministically spread by the FFT to 
neighboring subchannels. As a result, instead of one (or a few) subchannels being 
effected by the narrowband impairment, many subchannels are effected. The resulting 
loss of signal-to-noise ratio (SNR) can unacceptably limit the capacity of the 
communication channel 24. With continued reference now to Figure 2, the data is S/P 
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converted as shown by element 30, into blocks of length N+ P 9 subsequent to time 
domain equalization 26. The data is then transformed (e.g., in the case of ADSL, by 
removing the cyclic prefix as shown by element 32 and a FFT as shown by element 34) 
to length N. Following this transformation process, the present multichannel FEQ 12 is 
applied to each block of data. The multichannel FEQ 12 output data for the n th 
subchannel can then be expressed as 

x(n,t) = Y(k(nlt) T z(n). 

as discussed herein before with reference to equation (7). The output data from the 
multichannel FEQ 12 is the P/S converted, as shown by element 36, to form an estimate 
of the original data. 

Figure 3 is a flow diagram illustrating a method 100 of estimating multichannel 
FEQ coefficients for the n th subchannel, g(w), according to one embodiment of the present 
invention. As discussed above, the present method relates to T blocks of transmitted 
data, X(;, t) 9 t=l,...,T > md T blocks of received data, Y(:, t),t=l 9 ..., T depicted by 
elements 102 and 104 respectively. The received data can then be used to determine a set 
of subchannels, k(w), for the multichannel FEQ 12 to use in conjunction with subchannel 
n, for each subchannel n used to transmit data, as shown by element 106. The set of 
subchannels, k(w), includes the n th subchannel and one or more of the following: the n th 
subchannel's neighboring subchannels; subchannels where RFI is located; and 
subchannels where the noise power is (relatively) strong. Finally, the multichannel FEQ 
coefficients for the n ih subchannel, g(«), are generated as shown by element 108, via 
equation (2) using all of the received data at once, according to equation (9), or 
adaptively, using a block of data at a time, according to equations (1 1) and (12). These 
multichannel FEQ coefficients are then used in a complex multiply operation to adjust for 
the gain and phase of the effective channel response in each subsymbol. 

In view of the above, it can be seen the present invention presents a significant 
advancement in the art of communication between subscribers over shared media. 
Further, this invention has been described in considerable detail in order to provide those 
skilled in the data communication art with the information needed to apply the novel 



Patent Application of Arthur J. Redfern 



Page 12 of 18 



Docket No.: TI-3 1367 



principles and to construct and use such specialized components as are required. In view 
of the foregoing descriptions, it should be apparent that the present invention represents a 
significant departure from the prior art in construction and operation. However, while 
particular embodiments of the present invention have been described herein in detail, it is 
to be understood that various alterations, modifications and substitutions can be made 
therein without departing in any way from the spirit and scope of the present invention, as 
defined in the claims which follow. 
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WHAT IS CLAIMED IS: 



L A method of canceling communication system noise interference, the method 
comprising the steps of: 

a) receiving J blocks of data, Y(:, t) 9 t = 1, . . T, comprising T blocks of data, 
X(:,0, t = 1, . . T 9 transmitted over predetermined subchannels; 

b) determining a set of subchannels, k(«), for the multichannel frequency domain 
equalizer (FEQ) for subchannel n; 

c) generating multichannel FEQ coefficients, g(w), for the n th subchannel used to 
transmit the data; and 

d) performing multichannel (FEQ) for subchannel n using the generated 
multichannel FEQ coefficients. 

2. The method of canceling communication system noise interference according to 
claim 1 wherein steps b-d are repeated for each subchannel n used to transmit the T 
blocks of data. 

3. The method of canceling communication system noise interference according to 
claim 1 wherein the step of determining a set of subchannels, k(«), for a subchannel n 
used to transmit the Tblocks of data includes selecting subchannel n. 

4. The method of canceling communication system noise interference according to 
claim 3 wherein the step of determining a set of subchannels, k(w), for a subchannel n 
used to transmit the T blocks of data further includes selecting neighboring subchannels 
to subchannel n. 

5. The method of canceling communication system noise interference according to 
claim 3 wherein the step of determining a set of subchannels, k(w), for a subchannel n 
used to transmit the T blocks of data further includes selecting subchannels where radio 
frequency interference is located. 
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6. The method of canceling communication system noise interference according to 
claim 3 wherein the step of determining a set of subchannels, k(«), for a subchannel n 
used to transmit the T blocks of data further includes selecting subchannels having 
predetermined noise characteristics. 

7. The method of canceling communication system noise interference according to 
claim 1 wherein the step of generating multichannel FEQ coefficients, g(n\ for 
subchannel n 9 comprises solving the equation g(n) = Y(ny l x(n) , where Y(n)~ l is the 
pseudoinverse of a matrix of received data for subchannels k(w), and x(«) is a vector of 
transmitted data for subchannel w. 

8. The method of canceling communication system noise interference according to 
claim 7 wherein g(n) is determined adaptively using a block of received data at a time 
according to an equation defined by: g(n) - g(w) + ^(t)e(t)Y(k(n%t)* , where g(n) is the 
vector of multichannel FEQ coefficients for subchannel n and Y(k(/i), 0* is the conjugate 
of a matrix of received data for subchannels k(n). 

9. The method of canceling communication system noise interference according to 
claim 8 wherein values ofe(t) are determined according to an equation defined by: 

e(t) = X(n 9 t) - Y(k(n),t) T g(w) , where X(n, t) is the transmitted data for subchannel n at 
time t, and Y(k(n),0 r is the transpose of a matrix of received data for subchannels k(«). 

1 0. The method of canceling communication system noise interference according to 
claim 8 wherein pi(t) controls the adaptation according to least mean squares and has a 
value determined according to an equation defined by: R = 

where Y(k(n),t) H is the conjugate transpose of a matrix of received data for subchannels 
k(n). 
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1 1 . The method of canceling noise interference according to claim 8 wherein ju(t) 
controls the adaptation according to normalized least mean squares and has a value 

determined according to an equation defined by: a(t) = , 

P + Y{Un\t) H Y{Hn\tY 

where a € (0,2) , 0 < {3 , and Y(k(n),£) is a matrix of received data for subchannels k(n). 

12. The method of canceling noise interference according to claim 8 wherein ju(t) 
controls the adaptation according to power normalized least mean squares and has a value 

a 



a determined according to an equation defined by: /u(t) = 
where a 2 (t) = ccr 2 (*-!) + \e(t)\ 2 , c € (0,1) , and 0 < a < 



<y\t) 

_2_ 
M ' 



13. A system for canceling communication system noise interference, the system 
comprising: 

a multichannel frequency domain equalizer configured to receive T blocks of data, 
Y(:, t),t= 1, T, comprising Tblocks of data, X(: 9 f) 9 t= 1, r, transmitted over 
predetermined subchannels, wherein the multichannel frequency domain equalizer is 
operational to generate multichannel frequency domain equalization (FEQ) coefficients, 
g(n), associated with the n th subchannel used to transmit the data, and to perform 
multichannel FEQ for the n th subchannel using the generated multichannel FEQ 
coefficients, and further wherein the FEQ coefficients are associated with a set of 
subchannels, k(«), for the n th subchannel used to transmit the T blocks of data. 



14. The system according to claim 13 wherein the FEQ is operational to increase a 
subchannel signal-to-noise ratio beyond that achievable using a single channel FEQ. 



15. The system according to claim 13 wherein the FEQ is operational to cancel 
correlated subchannel noise caused by deterministic noise spreading associated with a 
plurality of subchannels. 
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16. The system according to claim 13 wherein the FEQ is operational to increase 
achievable digital subscriber network data rate beyond that achievable using a single 
channel FEQ. 
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Abstract of the Disclosure 

A multichannel frequency domain equalizer improves the subchannel signal to 
noise ratio by canceling correlated noise that is caused by deterministic spreading at the 
receiver that introduces correlation into the subchannel noise vector. 
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Transmit T blocks of data: 

xc,/y=i,-,r 



J? 



Receive T blocks of data: 

YU),t = \,-,T 



!o4 



Determine a set of subchannels, k{n), for the multichannel FEQ 
to use in conjunction with subchannel n, for each subchannel n 
used to transmit data, wherein the set of subchannels includes: 
i. the n lh subchannel and one or more of the following: 
h\ the n th subchannel's neighboring subchannels; 

iii. subchannels where RFI is located; and 

iv. subchannels where the noise power is (relatively) strong. 



Solve for the multichannel FEQ coefficients for the 
subchannel, g(n), according to: 

G(n,:) = argmin£b(«,0-G(w,:)Y(:,0| 2 J 

CO.;) " 

using all of the data at once, according to: 

g(H) = ¥(»)-' X(»). 

or adaptively using a block of data at a time, according to: 

g(«) = gO) + MfM0Y(k(«),f)* , where 

e(0 = ^(",0 _ Y(k(«),O r g(n). and where 
jU(t} controls the adaptation. 
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